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Abstract 

With the augmentation of IceCube by DeepCore, the prospect for detecting dark matter anni- 
hilation in the Sun is much improved. To complement this experimental development, we provide 
a thorough template analysis of the particle physics issues that are necessary to precisely inter- 
pret the data. Our study is about nitty-gritty and is intended as a framework for detailed work 
on a variety of dark matter candidates. To accurately predict the source neutrino spectrum, we 
account for spin correlations of the final state particles and the helicity-dependence of their decays, 
and absorption effects at production. We fully treat the propagation of neutrinos through the 
Sun, including neutrino oscillations, energy losses and tau regeneration. We simulate the survival 
probability of muons produced in the Earth by using the Muon Monte Carlo program, reproduce 
the published IceCube effective area, and update the parameters in the differential equation that 
approximates muon energy losses. To evaluate the zenith-angle dependent atmospheric background 
event rate, we track the Sun and determine the time it spends at each zenith angle. Throughout, 
we employ neutralino dark matter as our example. 



I. INTRODUCTION 



Astrophysical data require about 20% of the energy density of our universe to be in the 
form of unseen matter. The nature of this dark matter (DM) is a mystery that is key to 
resolving several problems in astrophysics and cosmology. A possibility is that DM consists 
of stable or very long-lived Weakly Interacting Massive Particles (WIMPs). However, there 
is no such particle in the spectrum of the Standard Model (SM) of particle physics and a wide 
range of new physics models have been proposed to introduce WIMP candidates. The most 
popular scenarios include supersymmetry (SUSY), extra dimensions, and n-plet extended 
models {n > 1). 

SUSY with 7^-parity [T] is a well motivated possibility that alleviates the hierarchy prob- 
lem, realizes gauge coupling unification and facilitates the seesaw mechanism of neutrino 
mass generation when augmented with right-handed neutrinos. 7^-parity conserving SUSY 
provides the lightest 7^-odd particle as a WIMP candidate. The cosmic microwave back- 
ground data from Wilkinson Microwave Anisotropy Probe (WMAP) [2] pins down the relic 
dark matter abundance to be QuMh"^ = 0.1123 ± 0.0035, which stringently constrains the 
SUSY parameter space. The popular minimal supergravity model [3] (mSUGRA) is defined 
by a set of only five parameters. 

Indirect searches look for cosmic ray excesses in the diffuse background or from point 
sources due to DM annihilation or decay into SM particles. Due to their high penetration 
ability, neutrinos help with detecting DM deep inside gravitational wells that include nearby 
sources like the Sun, the Earth and the galactic center j4H6|. DM-induced neutrinos from 
the Sun can be observed if the signal rate is competitive with the atmospheric neutrino flux, 
which is created by collisions of cosmic protons and nuclei in the atmosphere and is the 
dominant background below a TeV. 

Among the many experiments searching for high energy neutrinos, we focus on the Ice- 
Cube (IC) detector [7] that is capable of observing neutrinos with energies above 100 GeV. 
The installment of DeepCore [8] (DC) significantly lowers IceCube's energy threshold and 
enhances the ability of detecting neutrinos from light WIMP annihilation. We carry out a 
detailed simulation of IC/DC detection of the neutrino signal from neutralino annihilation in 
the Sun for the sample relic-density-consistent mSUGRA points listed in Table |I] in standard 
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TABLE I. Sample points in regions of mSUGRA parameter space that are compatible with the 
dark matter relic abundance; the sign of the /x parameter is positive and the top quark mass is 
172.7 GeV. The lightest neutralino is the WIMP candidate. Points A— D are selected to have 
m^o ~ 10^ GeV and large annihilation rates from parameter scans with tan /3 = 10 and 50. Points 
E— G are representative relic density compatible points from Ref. [£J. Masses are in GeV, and the 
last column gives the number of annihilations per year in the Sun. Only channels with branching 
fractions larger than 10"'^ are listed. 



notation; for a review of the various parameter regions see Ref. [TO] . 

In Section [IT] we summarize the physics of DM condensation and annihilation in the Sun. 
In Section |III| we describe our calculation of the spin correlated neutrino source spectrum 



and its propagation from the Solar core to the Earth. In Section IV we discuss the simulation 
of neutrino-induced events at the IC/DC detector. We present our results in Section [V| and 



summarize in Section VI In 5 appendices, we provide several details of our calculations. 



II. DM CAPTURE AND ANNIHILATION 

As the Sun sweeps through the dark matter halo, WIMPs collide with solar nuclei and 
become gravitationally trapped. The capture over a long period of time leads to condensation 
of low-speed WIMPs in the center of the Sun. The capture rate [5] Cc receives contributions 
from spin-independent (SI) and spin-dependent (SD) scattering between WIMPs and nuclei. 
Then, Cc = C^^ + C^P with [II] 
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SI / \ 

Cg^ = 4.8xlO^V^^^5^F./.0.^^(^) , (1) 
C^^ = 1.3 X 1029^-1^^4^^ ( ^1 , (2) 

where i sums over the elements with significant abundance in the Sun ranging from hydrogen 
to iron, po.3 is the local DM halo density in units of 0.3 GeV/cm^, V270 is the average DM 
dispersion velocity in units of 270 km/s, 171^. denotes the mass of the nucleus of the i^'^ 
element in GeV, and o"j is the SD/SI scattering cross section in pb. /j, Fi and 5* are the mass 
fraction, kinematic suppression and form-factor suppression |,12j for nucleus i, respectively. 
(pi describes the distributions of the i^^ element. We refer interested readers to Ref. [IT] 
for a detailed discussion and the values for these parameters. For most mSUGRA points 
consistent with the measured relic density, is greater than a^^ by two to three orders of 
magnitude, but does not necessarily dominate the capture rate. 

As the density builds up in the center of the Sun the annihilation of DM particles occurs 
more frequently. Eventually equilibrium sets in, Cc = 2Ca, where Ca_ is the annihilation 
rate. However, it was pointed out in Ref. [T^ that large areas of SUSY parameter space do 
not saturate this equilibrium condition and Ca can be significantly below Ccf^. The DM 
annihilation rate in the Sun [Tl] and can be parametrized by [llj 

Ca = ^tan^ (t/r) , 

i (3) 



t/T = 330 



s-1 cm3s-i l 10 GeV 



where t and r denote the age of the Sun and the equilibrium time scale, and {aAv) is the 
annihilation cross section averaged over the velocity distribution in the nonrelativistic limit. 
We do not assume that equilibrium holds, and calculate the annihilation rate for each region. 

III. NEUTRINO SOURCE SPECTRA AND PROPAGATION 

The source neutrino/antineutrino flux is 

^ = C^yBF.M, (4) 
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where i denotes each annihilation channel, and BFj and are the corresponding branching 
fractions and normalized (to each annihilation event) neutrino energy spectra, respectively. 

The dominant annihilation channels for our sample mSUGRA points are provided in Ta- 
ble [ij Note that the evolution of the mSUGRA renormalization group equations (RGBs) to 
the weak scale can be numerically sensitive to the GUT-scale parameters and lead to signif- 
icant uncertainty in the annihilation rate. This is especially true for the Focus Point region. 
We use DarkSusy [I5] to calculate {cav) and Cp^n^^, needed to determine the annihilation 
rate. 

Since neutralinos are Majorana fermions, their annihilation into light fermion pairs is 
helicity-suppressed in the nonrelativistic limit. Thus WIMP-induced neutrinos do not have 
a line spectrum in mSUGRA. The dominant neutrino source is the annihilation into gauge 
bosons, 2""^ and 3'''^ generation fermions, and their subsequent decays. Spin-correlation effects 
are visible in the neutrino energy spectrum, especially in the case of gauge bosons for which 
final states with transverse polarizations dominate in the static limit. It is important to carry 
out a spin-correlated calculation since the transversely polarized WW channel produces a 
significantly harder neutrino spectrum than the longitudinal WW channel. Similarly, a left- 
handed r produces more neutrinos than a right-handed r. In our analysis, we retain the spins 
of particles that directly result from the decay of particles pair-produced in DM annihilation. 
Secondary neutrinos arise from subsequent decays, and we include the spin correlation in 
helicity-dependent r decays. At this level the spin information of the primary and major 
secondary neutrino contributions are taken into account. The s-channel top-pair final state 
can be significant for large neutralino masses, and the polarization of the on-shell W from 
top decay affects the leading neutrino distribution. For this particular channel we proceed to 
the next level and keep the spins of the decay products of the W. See Appendices [A] and [B] 
for a discussion of our simulation of the various annihilation channels. 

The dense solar matter absorbs all the muons and relatively long-lived hadrons. A large 
fraction of c, b hadrons also scatter before they can decay. Thus muons are considered the 
end of the cascade and do not contribute to the neutrino flux. The absorption probability 
of c, b hadrons is determined by the competition between scattering and decay rates, which 
is discussed in Appendix [Cj 

The low flux density at the center of the Sun means neutrinos are created incoherently in 
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the flavor basis. Following Ref . [16] , we treat the neutrino propagation through the Sun with 
the flavor-density matrix that denotes the distribution in the flavor basis = Vei ^fi, ^r)- 
The propagation is governed by the equation, 



— = -z [H pi - — 
dr ' dr 



(5) 

NC^CC 



where r is the distance from the center of the Sun. The 4^ term denotes the neutrino flux 

dr 

attenuation from neutral-current (NC) and charged-current (CC) scatterings off the solar 
matter, and 're-injection' at lower energy after NC scattering, as well as secondary z/g, Vfj, 
production from r regeneration; see Appendix |D] for details. As neutrinos are created in 
gauge eigenstates, p is a diagonal matrix at the center of the Sun with diagonal elements pu 
denoting the fractions of the corresponding flavor. The flavor-basis Hamiltonian includes a 
rotation from a diagonal mass basis matrix and a matter effect term due to CC scattering 
with electrons [T7], 

H = -^V diag (0, 5ml^, ^m^,) ± v^G^nediag(l, 0, 0) , (6) 

where is the electron number density, Gp is the Fermi constant, and Smfj = mf — m| are 
the neutrino mass-squared differences. The sign of the matter term is positive for neutrinos 
and negative for antineutrinos. V is the neutrino mixing matrix parametrized by three 
mixing angles 6ij and a CP phase. We set the oscillation parameters to be |18j : 

Sml-^ = 8.1 X 10"^ eV2 , 5ml^ = 2.2 x lO'^ eV^ , = 33.2° , ^13 = , ^23 = 45° . 

(7) 

After leaving the surface of the Sun, neutrino propagation is dictated by the vacuum 
Hamiltonian ^q. Our choice 613 = causes the very long wavelength modes to be 
suppressed so that IceCube measures the average neutrino flux in half a year. The vacuum- 
oscillation average is obtained by dropping the off-diagonal terms of Vp, i.e., the neutrino 
density matrix in the mass basis. We ignore the attenuation of neutrinos as they pass 
through the Earth. Due to the low density of earth matter, attenuation becomes signiflcant 
only above 10^ GeV and is negligible for the atmospheric background and sub-TeV DM 
neutrino sources. Fig. [T] shows propagation effects on the neutrino spectra for Point A. 
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FIG. 1. Neutrino spectra for Point A at the center of the Sun (dashed) and at the detector (sohd). 
The dominantly transverse WW channel leads to a hard spectrum. The attenuation is mainly 
due to CC scattering in the solar medium. NC scattering and r-regeneration feed neutrinos back 
to lower energy. The spectra at the detector are time-averaged by removing vacuum oscillatory 
components. The and v-t spectra at the detector are almost identical. 

IV. DETECTOR SIMULATION 

The IceCube detector is a km^ sensor array that tracks muons. By selecting events that 
come from below the horizon, the only source of muons are the atmospheric and cosmic 
neutrinos that penetrate the bulk of the Earth and CC scatter with nuclei inside or in the 
vicinity of the detector. The detected muons are grouped into two categories: 'contained' 
muons with tracks starting within the instrumented volume, and 'up-going' muons that are 
created under the detector and range into the detector |Tn]. Deep Core is an extension with 
six additional strings inside the IceCube array. DC has a muon detection threshold as low as 
10 GeV and vetoes all the muons detected by the surrounding IceCube strings, thus elimi- 
nating most of the huge downward muon background and allowing Air detection of contained 
muons events. Recently Ref. [20] pointed out that cascade events may be detectable and 
may increase the event count at IceCube significantly. However, the angular resolution of 
cascade events is much poorer than track-like muon events, leading to a significantly larger 
acceptance cone size for atmospheric neutrinos. We found that the atmospheric cascade 
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event rate completely overwhelms the solar DM signal by a factor of ~ 10^ even with an 
optimistic angular cone size of 30°. Consequently, we limit our study to track-like events. 

The muon threshold energy has a significant effect on the background rates. Although 
IC and DC can detect muons with energy as low as 50 and 10 GeV respectively, the angular 
reconstruction for the muon track requires triggering of least three optical modules which 
raises the energy threshold. In most of what follows, we assume the threshold energy to be 
100 GeV and 35 GeV for IC and DC, respectively, and the angular resolution to be 1° half 
apex angle. 

The contained muons are detected at their initial (maximum) energy and the rate is given 

by 

— - ViEMO.) dE^ n„/, — — , (8) 

where is the differential cross section of creating a muon of energy E'^ from CC scattering 
off a neutron/proton, nn/p is the numerical density of neutrons/protons in the medium and 
r]{6z) is the detection efficiency at zenith angle 6z- Note that at the South Pole the Sun stays 
within the range of ~ 23° from the horizon. We optimistically assume that the efficiency 
has a weak angular dependence and set r/ = 1. This gives the maximum muon count which 
can be further adjusted with realistic detector information. is the incoming neutrino flux 
at the detector. V is the detector volume which we take to be an energy-independent 1 km^ 
for IC. For DC we parametrize the effective volume at SMT3-trigger level (DC veto) [21] as 



Vr 



DC 



0.32 + 7Mx + 8.91a;2 , < x ^ 1.09 

11.9 + 7. 77x- 1.06x2 , 1.09 <x ^3.0 (9) 
26.1 , X > 3.0 



where x = \ogiQ{E^/GeV) and Vdc is in megatons of water. The effective volume in ice is 

~ ■ Pwater/ P ice- 

The up-going muons lose energy before reaching the detector modules. Muons with energy 
below a TeV lose energy mainly via ionization. Above a TeV, radiative losses becomes 
significant. Following Ref [22], the up-going muon rate is given by. 



dE^dn 



A^{E„e,) / dz / dElP{E%E^-z) 
Jo Jeu 



dE^ Y: n^lv J^, (10) 

where E^^ and E^J^ denote the muon energy before and after propagating a distance z outside 
the detector. P{El,E^ ; 2) is the survival probabihty of the muon after propagation which we 
simulate using the Muon Monte Carlo package [22]; for details see Appendix [Ej A^^E^.O^) 
is the effective muon detection area (22] parametrized by 

A^ = Ikm^ ■ Ao{E^) (0.92 - 0.45 cos 9,) (11) 

I , X < 1.6 

Aq = I 0.784(x - 1.6) , 1.6 ^ a; < 2.8 
0.9 + 0.54(x-2.8) , X ^ 2.8 

where x = logio(-E'^/GeV). 

Observation of up-going events at IC starts with the September equinox [t = 0) and ends 
with the March equinox {t = 0.5). During this interval, the zenith angle of the Sun follows 

^^(t) = 90° + 23.43°sin(27rt) (0 ^ t ^ 0.5) . (12) 



For IC, the muon rates are obtained by integrating Eq. [8] and Eq. 10 over the real-time 
zenith angle of the Sun; DC operates year-round so that observation time doubles. While 
the solar core is a point-like source, the major background for the DM signal comes from 
atmospheric neutrinos, which depends on the angular resolution of the detector. The direc- 
tional atmospheric neutrino flux is measured by Super-K [23] and is symmetric about the 
horizon (6z = 90°). The left panel of Fig. |2] assumes a cone of 1° (half the apex angle) and 
shows the muon rate from atmospheric neutrinos at IceCube and DeepCore as a function of 
zenith angle. From the right panel it should be noted that the background rate increases 
quadratically with the opening angle of the acceptance cone. For a more detailed description 
of up-going and contained atmospheric background events, see Appendix C of Ref. 



V. DARK MATTER SIGNALS 



The number of energetic neutrinos above detector threshold determine the prospects for 
detecting new physics at IceCube/DeepCore. For a solar WIMP signal there are three major 
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FIG. 2. Left panel: Integrated atmospheric background rate as the Sun sweeps through each degree 
in zenith angle in half-a-year at IceCube and a full-year at DeepCore assuming a 1° acceptance 
cone. Right panel: The acceptance cone-size 69 (half the apex angle) dependence of the yearly 
atmospheric background rates. 



10 



24 



23 



10 



^ 10^' 

G 
G 

< 

21 

^ 10 



10 



20 



I e 
• • 



OA„..B 



50 100 150 200 250 
[GeV] 



300 



FIG. 3. Annihilation rate for mSUGRA points that are consistent with the measured relic density 
at the 95% C. L. Plotted points are from two scans in the (rriQ, mi/2) plane with tan/3=10 and 
50. The light gray dots represent parameter points excluded by XENONlOO/Super-K/IC data. 
Points A-D of Table |l] are marked among the nonexcluded points with tan /3=10 (blue circles) and 
tan/3=50 (red triangles). 
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TABLE II. Atmospheric background rate and signal rates for the points of Table |Ij The observation 
time for IC and DC are ^ and 1 year, respectively. The acceptance cone has a 1° opening angle 
(half the apex angle) for both background and signal. Ej^j^^, denotes the muon energy threshold. 

factors: (z) annihilation rate; (ii) muon energy threshold vs. WIMP mass; {Hi) annihilation 
channels that produce energetic neutrinos. 

The relic density provides a stringent constraint that relates the first two factors. Fig. |3] 
illustrates the dependence of the yearly annihilation rate on the neutralino mass in different 
regions of mSUGRA parameter space. The inverse dependence of the annihilation rate on 
WIMP mass is evident with a minimum mass set by the energy threshold of the detector. 
Viable mSUGRA points must also be consistent with the XENONIOO [26] constraint on asi, 
which is the most stringent among nuclear recoil experiments |27] , as well as the Super-K |28j 
and IC [23 constraints on asn- The lower energy threshold of DeepCore greatly enhances 
the signal with respect to the background and can be crucial for WIMPs with mass ~10^ 
GeV. The Focus Point region allows a large aso coupling and is the most popular discovery 
scenario for mSUGRA. The bb channel generally produces less and softer neutrinos compared 
to the WW and ZZ channels. The t~^t~ channel provides energetic neutrinos but often has 
a low branching fraction. The neutrino signal from a massive WIMP becomes hard to detect 
as the bb channel dominates (Point D). 

In Table [Tl| we list the signal rate for the points in Table jl], and the atmospheric back- 
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FIG. 4. Time-integrated muon energy spectra for the atmospheric background and WIMP signals. 

ground rates; for IC we have entertained two possible detector thresholds for comparison. 
The differential energy spectra for Points A and D are shown in Fig. |4j With sufficient 
statistics as for Point A, it is even possible to construct the shape of the energy spectrum. 



VI. SUMMARY 

We presented a calculation of neutrino source spectra from solar dark matter annihila- 
tion, accounting for spin-correlations, and a simulation of the propagation and detection of 
these neutrinos at the IceCube/DeepCore detector. We considered the mSUGRA model to 
illustrate spin-correlations, but similar techniques can be applied to other WIMP models. 

With an angular resolution of 1° half apex angle and muon energy threshold of 100 GeV 
and 35 GeV for IceCube and DeepCore respectively, the yearly atmospheric background rate 
is 5.6 for IC upgoing events, 14 for IC contained events and 2.3 for DC contained events. 
Generically, a neutralino annihilation rate of 10^^ yr~^ is necessary to compete with the 
atmospheric background. With its lower energy threshold DC plays an important role in 
detecting neutrinos from relatively light (~ 10^ GeV) dark matter candidates, which have 
less suppressed annihilation rates. For DM masses that are not much above the detector 
threshold, accounting for the helicity distribution of the final state particles can be critical for 
the detectability of the signal. For example, W pairs produced from neutralino annihilation 
are transversely polarized and give a much harder neutrino spectrum than if both helicities 
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contribute equally. 

The neutrino-copious W~^W~ ,t^t^ and ti channels are the main contributors to the 
neutrino signal. Since neutralinos are Majorana fermions, spin correlation requires W^W~ 
to be transversely polarized which yields a hard neutrino spectrum and enhances the muon 
rate above IC/DC thresholds. The Focus Point region is the primary mSUGRA discovery 
region for IC/DC because of the high annihilation rate. In general, the neutrino signal rate 
can be compromised if the annihilation occurs primarily into channels that do not produce 
a significant neutrino flux or a hard spectrum. 
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Appendix A: Source neutrino spectrum 



The primary channels that contribute neutrinos are given in Table III Pre-shower spectra 
are generated with MadGraph/MadEvent (MG/ME) |30j or Calchep [3l| that keep the spin- 
correlation of the final state particles. We modified the MG/ME-Pythia [321 133] interface 
to develop the shower and take into account the solar absorption of 6, c hadrons. Secondary 
neutrinos come from subsequent decays of taus and 6, c hadrons in the shower. In the case 
of r decay, helicity information is facilitated by the Tauola [31] package, as part of the 
MG/ME-Pythia interface. 

In the W^W~ and ti channels, the undecayed W and t are made stable in the shower 
and their contribution is added as the charge conjugate of the neutrino spectra. Similarly, 
in the ZZ channel one Z boson is tagged stable in the shower and the neutrino spectrum is 
doubled. The r^r^ channel needs special treatment as a limitation of the LHE format ^35j 
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Channel 


Final state 


w+w- 


W+*UrT-, W+*sc, W+*dc 


ti 


t*br>rT~ , t*bsc, t*bdc 


ZZ 


Z*T+T-, Z*u^u7, Z*bb, Z*cc 


Zh 


h*T^T^, h*Vj^Vj^, h*bb, h*cc, 




Z*T+T-, Z*bb, Z*cc 






bb 


bb 


cc 


cc 



TABLE III. Neutralino annihilation channels that contribute neutrinos. Particles marled with a * 
are made stable in the shower. The spin correlation in the r^r^ channel is treated separately with 
helicity-dependent decays. 

makes it difficult to pass the s/t/u channel particle information to Pythia, and Tauola cannot 
reconstruct the helicity of the pair-produced r leptons. We circumvented this problem by 
treating r decay analytically for the t^t~ channel, as discussed in Appendix [b] 

MadGraph has difficulty in producing transverse-W^ spectra for a few mSUGRA points 
in the static limit and we switched to Calchep for the W^W~ channel. A drawback is that 
Calchep sums over final spins, so helicity information of the r leptons is lost. Thus, most 
of the secondary (soft) component in the W^W~ channel is obtained from unpolarized 
r decays. In the ti channel, final state radiation is turned off in Pythia to avoid a problem 
of flavor sum of parton clusters, after t is tagged stable. The resultant h quark energy may 
be high by a few percent. However, the effect is irrelevant since secondary neutrinos are 
dominantly produced from r decay. 



Appendix B: r decay 

While subdominant when the W~^W~ , ZZ and tt channels are kinematically allowed, the 
r+r^ channel is a major source of neutrinos for lighter DM as the other dominant channel, 
bb, produces less neutrinos per annihilation. The relevant decay channels are r — )■ Ut-Wi and 



1^-,-+ hadrons as listed in Table IV The neutrino energy spectrum in the rest frame of the r 
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can be expanded as 

where x = 2Ey/m^ is the energy fraction, and /o, /i are projections of the distribution on 
the first two spherical harmonics. After boosting into the lab frame the first two harmonic 
coefficients are 



9Q{y) = / dxfo{x)/x, 



gM= / dx{2y-x)h{x)/x\ (B2) 



where y = E,^/Et-. The lab frame neutrino spectrum is 

1 dN 



No dy 



9oiy) + P9iiy), (B3) 



where P = ±1 for a left /right-handed r. For a two-body decay /o_i = S{1 — x — m\/m1) 
while for the r — )■ Vr^vi channel, these functions are 



2x^(3 - 2x) for z/^ -2x'^i2x - 1) for 

/o = <( , /i = , , (B4) 

12x (1 — x) for vi 12x (1 — x) for Ui 



and the lab-frame spectra are listed in Table IV For the short-lived mesons p and oi we 



smear the Dirac-5 to account for their widths. The rest-frame distributions are modified 
with a Breit-Wigner approximation, 

/o, = 5{l-x- r_) -> = C \ , (B5) 

V-L I mesj ~ '"'mes'- mes"''T 

where r^es = ^%ai/™Ty ^mes is the meson decay width and C is a normalization factor. The 
neutrino spectra resulting from left and right-handed r decays are shown in Fig. |5] 



Appendix C: Hadron absorption 

Hadrons that contain c, b quarks contribute to the neutrino fiux through their weak decay 
modes. However, the dense environment at the center of the Sun shortens the mean free 
path of these hadrons and absorption effects become significant when the mean free path 
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Decay mode (i^r) 


BF 




gi{y) 




0.18 


5 
3 


- 3l/2 + 4y3 


1 

3 


+ - 3y2 




UrTT 


0.12 


1 


9(1 -r^-y) 


(l-r,)2 - 


-y) 




0.13 


/' 






-2x)/i*(x)x- 


'^dx 


UrP 


0.26 


I' 






- 2x)/f 


^dx 




0.13 


Jy 


fx{x)x''^dx 









Decay mode {i>i) 


BF 


go{y) 


gi{y) 




0.18 


2 


- + iy^ 


-2 + 12y - 18y2 + 8y^ 



TABLE IV. T decay modes, their branching fractions and fragmentation functions go and gi. The 



smeared distributions /* are as in Eq. B5, We approximated fx{y) for the inclusive r — )• 
channel with a generic phase-space decay into 4 pions and Vt-. 




FIG. 5. i^Tji^i distribution from the decay of left (L) /right (R) handed taus. 



is comparable to the decay length. We take the nucleoli scattering cross sections of c, b 
hadrons [15j to be 

1.4 X 10^^° for mesons 
2.4 X 10"^° m^ for baryons 



a{E) 



(CI) 



and the mean free path is 



nQo{E) 



(C2) 



16 



cc channel 




'Irri n 

60 80 




20 40 60 80 
Ey [GeV] 



FIG. 6. Neutrino (solid) and antineutrino (dashed) spectra from the cc (left panel) and bb (right 
panel) channels for a 150 GeV WIMP. 

where Uq is nucleon number density at the center of the Sun. The decay length is 



XdiE) = cTj{E) 



(C3) 



where 7 is the Lorentz boost factor, c is Ught speed and r is rest-frame Ufe time. 

We assume that each scatter with a nucleus leaves the hadron with an average fraction 



0.7 for b hadrons, and e = 0.65- 



for c hadrons 



of the hadron's initial kinetic 



energy Eq. After the n*^ scattering, the ratio of the scattering probability to the decay 
probability is 



scattering 
Pdecay 



Xd{En) 
As (-En) 



(C4) 



where En ~ 



Eoe"'. We modified the Pythia package to simulate the decay and energy loss 
of c, b hadrons. The probability of scattering and decay are calculated iteratively for each 
hadron on an event by event basis. When a hadron is considered decayed, it is handed back 
to Pythia which develops the decay and showering. Since not all scatterings are elastic, we 
allowed the hadrons to scatter only once before decaying or being stopped/destroyed. The 
resultant neutrino spectra for a 150 GeV WIMP are shown in Fig. [6] 
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Appendix D: Neutrino scattering in propagation 



The neutrino flux becomes attenuated by scattering off solar nuclei as it traverses the 
Sun. NC scattering decreases the neutrino energy while CC scattering converts the neutrino 
into the corresponding lepton whose decay introduces a secondary influx of less energetic 
neutrinos. The latter is often called "regeneration" . While electrons and muons are quickly 
absorbed by the dense solar medium, r leptons decay promptly and contribute to the soft 
component of the neutrino flux. 

The flux attenuation caused by NC scattering is flavor blind, 

= - r^,/„</^(i?.)p(i?.), (Dl) 
where p is the flavor density matrix and np/„ is the solar proton/neutron number density 



The injection of scattered neutrinos is given by 



= np/n / dE^ — p(^J, (D2) 

NC J "-^t/ 



dr 



where the dummy E'^ integrates over energy above E^. 
The attenuation via CC scattering is given by 



dp{E,) {Tccp} 



dr 

where Tec = diag (np/„ ^^'^(E^)) 



(D3) 



cc ^ 

CCi 



and the secondary contribution from r regeneration is 



dpij{E^,) 



dr 



reg. 



CC 



= rip/ 



dE' 



dEr 

^CC,UE,Uf_i 
p/n 



dEr 

-) 

-g,^{E,/E,)p^,{El] 



da':'rf^''''{El,Er] 



dE^ 



g,^jE,/Er)prAEl] 



(D4) 



where 11; = SuSji is a 3 x 3 matrix with only the / diagonal element nonzero and / = 1, 2, 3 
for e, fi, T respectively. 11/ picks out the diagonal terms from p. p denotes the density matrix 
for antineutrinos. Note that since neutrinos have definite helicity, the r leptons from CC 
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scattering are polarized; r~ is left-handed while is right-handed. See Appendix |b] for the 
helicity fragmentation functions QuXy) for r decay. The neutrino-nucleon scattering cross 
sections are evaluated with the numerical package Nusigma [5^ . 



Appendix E: Muon propagation in ice 

The muon energy spectrum z) softens after propagation in ice, as generically de- 

scribed by 

/■^" da(E^,E') 
V, d^ME,, ^) = - / dE'^ME,, zHz)^^^ 

rEmax daiE" E ) 

+ / dE;<P,{Elz)n{z) +dEME,)ME,.t)), (El) 

where f ^ ~ c is the muon speed, a is the muon scattering cross section, n{z) is the target 
density, and a describes ionization energy loses. E^j_ is the muon's energy at a propagated 
distance z and dummy variables E'^, E'^ denote the energy below/above E^^. In terms of the 
survival probability P{E^, E^; z), the final energy spectrum can be written as 

MEf., z) = j P{E'^, E„ z - zo)ME',, z,)dE'^ , (E2) 

where is the initial muon energy. For 2; 7^ 0, PiE'^, E^; z) can be below unity upon 
integration over i?^ due to muons being stopped before reaching z; P{E^^, E^; z) can be 
obtained from Monte Carlo simulations. 

The spatial integral in Eq. [lOjcan be done separately to yield the so-called "effective muon 
range" |22], 

/•oo 

R{E%E,)= / P{E%E,-z)dz. (E3) 
Jo 

R{E^, Efj_) represents the average incremental distance a muon travels per unit energy loss 
(in GeV) at E^. 

The average neutrino effective area is [22] . 

A:ff{E^,9^) = - J2 jdE.dEln^,, ^' R{E%E,)A,{E,,e,) , (E4) 



where A^{E^,6z) is the muon effective area given by Eq. 11, The attenuation of sub-TeV 
neutrinos inside the Earth can be ignored. 

19 



5 

s ' 

3 

O 

S' 2 



1 



0^ 
0.0 




0.2 



0.4 



0.8 1.0 



0.1 



^ 0.01 



0.001 



10" 









IceCube 




MMC simulation - 




Ave. floss approx. 



200 



400 600 
Ey [GeV] 



800 



1000 



FIG. 7. Left pannel: Differential effective range R{E^,E^) in ice for 100 GeV (dashed) and 
1 TeV (solid) muons. Right panel: Comparison of Al^^ from the IceCube detector simula- 
tion [22j, our simulation with MMC, and the approximate method using a parameterization of 
average neutrino energy losses. The plotted Al/^ is an angular average over 9z > 90°. With 
a = 2.5 X 10^'^ GeVcm^/g, /3 = 4 x 10~^ cm^/g, the average E,^ loss method agrees well with 
sophisticated simulations. 

We used Muon Monte Carlo (MMC) to simulate the muon survival probability P{E^, z). 
The interpolated effective range for 100 GeV and 1 TeV muons are shown in Fig. [7j As a 
test we calculated A^J^ and compared to the full detector simulation, and found excellent 
agreement as is evident from Fig. [7} 

Alternatively, a frequently used approximation for muon propagation parameterizes the 
average muon energy loss, 

-p{a + PE^) , (E5) 



dz 

where p is the medium density and a, /3 account for ionization and radiative effects. This 
procedure ignores the smearing of the muon energy distribution during propagation. For 
ice, we found a = 2.5 x 10"^ GeV cm^/g, /3 = 4 x 10"^ cm^/g agrees well with the MMC 
results, as shown in Fig. [7| There is a mild degeneracy in the parameters a and /3, so that 
a = 3 X 10"^ GeV cm^/g, /3 = 3 x 10~^ cm^/g [25], works just as well. 
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